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^ ' Abstract. We present all Petrov type D pure radiation space-times, with or 

without cosmological constant, with a shear-free, non-diverging geodesic principal null 
congruence. 



a'. 

I . 1. Introduction 
5-1 . 

bJO, 

A space-time is said to belong to Kundt's class if it admits a null congruence, generated 

^ by a vectorfield k, which is non-diverging (i.e. which is both non-expanding and non- 

Q"^ ■ rotating and therefore also geodesic). In the case of a pure radiation field, which, by 

^ ■ definition, has an energy-momentum tensor of the form Tab = (pkakb, with kak°- = 0, 7^ 

' 0, one can furthermore show[2], making use of the null energy condition Rabk^k^ > 

^ ■ and the Raychaudhuri-equation 

o : . 



X 



that k and k are aligned and that the associated null congruence is shear-free : 



sj . Rabk^k" = -2aa ^ a = and Rabk^k" = 0, 

Rab = (pkakb + Agab =^ k (X k, 

where A is the cosmological constant. The Goldberg-Sachs theorem implies then that the 
Kundt space-times must be algebraically special (type II, D, III or N or conformally 
flat) and that k is aligned with a repeated principal null direction of the Weyl tensor. 

Conversely, for aligned Petrov type D pure radiation flelds, it follows immediately 
from the Bianchi equations (without invoking the energy conditions) and from 

^0 = ^1 = ^3 = ^4 = 0, ^2 ^ 0, (1) 

$00 = $11 = $01 = $12 = $02 = 0, $22 ^ 0, (2) 

that 

K = (J = A = 0, (3) 

such that the principal null direction k aligned with the radiative direction is geodesic 
and shear-free, i.e. the converse part of the Goldberg-Sachs is generalized to this situa- 
tion. Moreover, it is known that aligned Petrov type D pure radiation fields cannot be of 
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(null) Maxwell type [4] nor of (null) neutrino or scalar field type [3], while Wils[T] showed 
that k is necessarily non-twisting {p = p). The diverging solutions (p ^ 0) then belong 
to the Robinson-Trautman class and are all explicitly known [2]. In the non-diverging 
case p = 0, however, only a few type D examples are available[3]; here we present all 
solutions belonging to the latter class. 

As shown by Kundt[5] in 1961 the line-elements admitting a geodesic, shear-free and 
non-diverging null congruence can all be expressed in the form: 



in which P is a real function of (C, C? and H, W are respectively real and complex 
functions of {(,(:U,r), to be determined by appropriate field equations. The vacuum 
solutions of this type have been known for a very long time[6]. A procedure (based on 
Theorem 31.1 in p]) is also available allowing one to generate non- vacuum solutions 
of Kundt's class, from vacuum solutions. However, in this way the Petrov type of the 
metric generally will be changed from D to II. Insisting that the Petrov type does not 
change, constrains the function Ho of the 'background-metric' and it was not certain 
whether one could generate all pure radiation solutions in this way. This technique was 
used in [3], where the authors managed to construct a family of type D pure radiation 
fields and where they conjectured that these solutions were the only aligned type D 
pure radiation fields of Kundt's class. Below we show that the solutions obtained in [3] 
only cover a small part of the entire family, which we present in this paper. 

The most obvious way to find all type D pure radiation fields of Kundt's class, is 
to start from the general Kundt metric (Jlj), express that the solutions we are looking for 
are of Petrov type D, and to make sure that the field equations for pure radiation are 
satisfied. This however introduces a hard to solve system of non-linear conditions. We 
therefore prefer to start from the general Geroch Held Penrose (GHP)[7] or Newman 
Penrose (NP)[8] equations, where we follow the notation and conventions of p], and to 
extract from these all possible invariant information, before introducing any coordinates. 

2. GHP analysis 

Introducing the set of conditions ([I])-(l3l) and the non-diverging condition p = within 
the GHP Bianchi equations, one arrives at 



ds^ = 2P-MCdC - 2du(dr + Hdu + WdC + WdQ, 



(4) 



P^2 = 0, 

g^2 = 3r^2, 
= -37r^2, 

I>'^2 = -3p^2, 



P$22 = 0, 

a$22 = 317^2 + r^22, 
g'<l>22 = 3z/^2 + r$22. 



(5) 



The GHP Ricci equations then reduce to 
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and 



^^ = ^'' gV=(r-7r)z/, (6) 

o 

Pu = P'n + (r + 7r)/i, 

P'fi = du - fi'^ + Vn - UT - $22, (7) 

P>/i = Svr + TTTf + \&2 + — • 

3 

The remaining basic variables at this stage are /x, z/, vr, r, \l/2 and their complex 
conjugates, together with $22; the cosmological constant A appears as a parameter. 
From the integrability conditions for \l/2 one obtains : 



Pu = 0, = 0' 

p/i = TTTf - rr, = -{tt + "^2 + —), (8) 

3 

g/x = -PV, PV = -(r + 7r)/i. 

[9', f>] /i now yields the algebraic equation 

(^2rr + ^2 + ^^ 7r+ (^2rr + ^2 + ^^ r = 0. (9) 

From [9, we can extract an expression for Q'P't : 

a'pV = - 2^2 - 2rr - /i + (^^2 + 2rr + /I - ttPV. 

If one then calculates [5, S'] P'r, making use of [P',I>]r, [S, P']r, and ([9]) the following 
expression for P't is obtained, 

P't = JIt — jM — 2/ir, (10) 

after which [!>', 5] vr results in a second algebraic equation 

/xg = 0, g = 6(r7f + r7r) + 3(rr + 7r7f)-A + 3(^2+'*2). (11) 

From S(l9l) — 2r(l9]) one readily infers 

TTTf = rr. (12) 
In the case where vr = r = 0, Svr = and P'Svr = yield 

vl>2 = -|, /i = (13) 

and all the above equations are identically satisfied. Consider now the remaining case 
TIT 7^ 0, and introduce the zero-weighted (i.e. spin- and boost-invariant) quantities 

A = |7r7f|i/2 = I rr 1 = ^ > Q), (14) 

T = = = = {T = T-^). (15) 
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In this way vr and r are conveniently parametrized by 

71=^, r = ATB, (16) 
the weight of both being absorbed in the (l,-l)-weighted quantity 



Within this notation, the equations (|9]) and f[TT|l translate into 

(T + 1) (2^2 + A/3) + *2 + T^2 = 0, (18) 

fiQ = 0, Q = 2{T^ + T + l)A^ + T{^2+^2-A/S), (19) 

respectively. Furthermore, we define the (-2,-2)-weighted quantities 

X = du, N~Bv (20) 

and the zero-weighted derivative operators 

V = B-^^, V' = B&. (21) 

From (I5|)-(l7l), and [3,5']z/ one then deduces the following autonomous first- 
order system for the zero-weighted fields ^^2, ^2, T, A and the (-2,-2)-weighted fields 



X X, N, N, $5 



1)^2 = 3AT^2, ^^^2 = -3A^2, (22) 

VT = AT{T + 1), VA = -A^ - '^^1±A^ (23) 

6 

VX = A{AT + 1)X - ATX - 2A^T{T + 1)N + 2A^{T + 2)NA^ 

+ (3^2 + ^2 + A/3)iV-A(r+l)<l>22, (24) 

VX = A{T - 1)X + (4^2 + A/3 + 2A'^)N, (25) 

VN = X^^^^N, m=(^^-'^]N, (26) 



6A 

V^22 = AT^22 + 3^2iV. (27) 

For later use, one immediately observes from (l23ll that T constant implies T = — 1. 
Also notice that (l22l) . fl23l) forms a zero- weighted subsystem, whereas the equations 
(l24l) -(l27l) are (-2,-2)-weighted, such that their right hand sides are necessarily linear 
and homogeneous in X, X,A^, X, $22- The complex conjugates of (l22l)-(p7ll yield an 
analogous V system. However, as deduced from [S, 9'] $22, both systems are constrained 

by 

3T(^2X - ■^2X) + 6A(T2*2X - ^2X) + r(^2 - '^2)^22 = 0, (28) 

which is again linear and homogeneous in X, X, X, X, $22- 
We are now ready to prove: 
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Theorem W.r.t. a Weyl canonical tetrad ^),^ any aligned Petrov type D pure 
radiation field, for which the principal null direction, aligned with the radiative direction 
is non- diverging, satisfies the following boost- and rotation-invariant properties: 

K = a = p = 0, X = fi = 0, (29) 
7r + r = 0, = ^2, TT $22 + z/ |^^2 + = 0. (30) 

Proof It immediately follows from ((71) that p = u = is inconsistent with $22 7^ 0, 
and from ^ that vr + r = implies ^^2 to be real. By ( fT3l) the theorem is true in the 
case where vr = r = 0, and it remains to establish fi = and tt = — r (T = — 1 in the 
above notation) when vrr 7^ 0. 

Referring to (fTOl l. we first treat the case Q = 0. Two linear and homogeneous 
systems in X,X,N,N and $22 are obtained from [(EHD, V'l^, VV'I^] and 

PP(l28l). respectively VVIH^. As $22 7^ 0, the determinants of these systems must 
vanish. This yields two polynomial relations in A^, A, \E'25 ^2 and T; eliminating 
and A by means of ffTSi) and Q = 0, and appropriately scaling the results with ^^2, 
yields two polynomial relations in the zero-order variables T and = ^^2/^2 of the form 
(T + l)^T^Pj(T, k) = 0, i = 1, 2, where Pi and P2 are irreducible and non-proportional. 
Thus T is constant, whence T = — 1 and \E'2 = ^2- Substituting this into Q = yields 
TTTf — \E'2 + A/6 = 0, and calculating the P'-derivative hereof gives /i = 0. 

Thus /i = in any case, and (JT]) implies 

g'z7 = 5/7 = -7rI7- 7fz/ + $22, i.e., X = X = -A(iV - TA^) + $22. (31) 

Substituting this in [l^,Vl^,V'^], {^,V^,VVl^] and [(12811 .1)^(12811.^^(12811] 
yields three linear and homogeneous systems in N,N and $22- Eliminating and A 
from the three determinants and (fTSll leads to two polynomial relations in T and k with 
(T -|- 1)'^ as the only common factor, such that again T = — 1 and \E'2 = ^2- Inserting 
this and ( l3Tl l into ( l28l l implies N = N (i.e., u/n is real), a final P-derivative of which 
yields 

A$22 = Ar(^2 + A/3), 
the translation of the last equation of (l30l l. □ 



One checks that, given the specifications ([II), ((21) and (l29ll . (l30ll . the remaining GHP 
equations are consistent, such that corresponding solutions exist. Referring to ([TOll and 



—2/3 

the above proof, we finally remark that all derivatives of Q \E'2 ^^e zero, such that 
one has the following important relation: 

A 2/3 

TTTT = \I'2 — — + ci\E'2 , ci constant. (32) 

Because of ([T3ll . this only yields new information when vr 7^ 0; the case ci = 4^ Q = 
will be distinguished from Ci 7^ for the corresponding integration in § 13.21 
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3. NP analysis 

In order to construct the metrics, we now switch to the NP-formalism, where we 
first introduce the invariant information, as obtained above. Appropriate boosts and 
rotations allow one to put e = and (3 = — (a + vr). This fixes the null tetrad up 
to rotations and boosts satisfying D^^ = and DA = = 6A = 6A. As z/ 7^ and 
Dz/ = (P — 3e — e)z/ = 0, one can completely eliminate the rotational degree of freedom 
by requiring V = v. By ( l30l ) this implies vf = tt. The solutions can then be divided into 
two families : one in which vr = and one in which tt 7^ 0. We also introduce 70 and 
7i/z/, the real and the imaginary part of 7, respectively. 



3.1. The family vr = 

In this case (l8|) shows \E'2 = — |^ 7^ 0. The NP Ricci equations can then be rewritten in 
the form : 

$22 = -4az/, (33) 



and 



Da = 0, Di/ = 0, 

D7o = -^, (34) ~5v = -2au, (35) 

D71 = 0, -'2aiy, 

Slo = 0, (36) 
uAa - ihi = 0, (37) 

Sa + 6a = Aa'^ + ^. (38) 

The Bianchi-equations now yield an expression for 6a (by which (i38l ) becomes an 
identity) : 

4 

By (l33ll and (l35l) u cannot be constant and this allows one to use u as a coordinate. 
As is clear from (l34l ). 70 provides a second independent function and hence can be 
used to define a coordinate r = — 270/A. As ^70 = it follows that o;^ = dr + Huj^ 
{H = — Ar, a;* representing the basis one-forms). From (l35|) we see that we can write 
(jj^ + CL>^ as : 

uj^ + uj^ = - . 39 

2 av 

The Cartan-equations read 

da;^ = 2a LJ^ A + 2i7i/z/ cj^ A a;^ (40) 
du;^ = 0, 

du'^ = u{u^ + a;^) A u;^ - 270 a;^ A a;^. 
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showing that a function u exists, such that cj^ = du. A fourth coordinate x is defined 
by writing uj^ — u:^ = i^Ndu + Udu + 4z/Xdx) , with A^, U, X real functions of (m, r, x, z/) , 
where the factor u has been introduced for later convenience. From (l40ll it then follows 
that the functions N, U and X are independent of r, so a transformation of x can be 
used to put = 0. Hence we obtain the following basis one-forms : 

= du - - — dv + 2wXdx, 41 

= du, 
LJ^ = dr + Hdu, 

where U and X are real functions of (m, u, x) and if is a real function of (m, r, z/, x). The 
corresponding directional derivatives are then given by : 

d i d 
S = -2av— - 



du AuX dx ' 

9 rrd , d U d 
A = H— + Au— 



du dr du AuX dx ' 
dr 

Note that, because of the remaining boost freedom, r is fixed only upto transformations 
of the form r gi{u)r + g2{u). 

From the directional derivatives of a we see that a can be written in the form 



a/— 2Az/^ + sm[u 



2 



a = s = ±1. 42 

Au ^ ^ 

Applying [A, D] and [5, A] to r, we obtain an expression for H{u,r,u,x) : 
H{u, r, z/, x) = — + -J- + n{u). 

From [6, 6]x we get X{u, u, x) = X{u, x), while [6, A]x shows that 
dX 1 fdU 



du Au \ dx ' ^^'^'^ 



and 

dU U 7i 



du u au"^ 
From f l37ll and D71, we see that 71 = 7i(m, x) and 



(44) 



du Xsu \ dx 
Making use of the latter equation, we can write [5, A]z/ as 
dAu/u 



sm^XAu - Aau^ . (45) 



dx 

and 

dAu/u d'^ijdx 
du ~ ~ Aau^X 



-167iaX (46) 

(47) 
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We now distinguish two subclasses m = or m 7^ (in both classes h and n are 
arbitrary real functions of u) : 

3.1.1. The class m{u) = 

Notice that, if m{u) = 0, ([42]) shows that A < and we put ao = a = y/-2A/A. ^ 
shows that 71 is independent of x. (l44ll . (I43l) and (l47l) then yield 

/\u{u,u,x) = (48) 



and 



[/(u,z/,x) = -^ + z/f/i, (49) 
2aQU 

where Ui is a function of (m, x). 
As a; ^ is now given by 

= -\^^^^^ + wUi\du ^—du + 2iuXdx, (50) 

2 \ 2aQi' J Aa^v 

with 7i a function of u and f/i,X functions of (u,x), we see that ii/f/i(M, x) can be 
absorbed by a transformation of x. The remaining equation to be satisfied fl46l) then 
reads 

d'^X dXdX . , , 

-—X - — — - 167i«o^' = 0, 51 
axau ou ox 

where there is still freedom in the choice of x. 

If 7i = the solution of (ISTil is of the form X(m,x) = f{x)/b{u), where we can put 
/(x) equal to 1 by a transformation of x. After a transformation of the r, z/ and u co- 
ordinates, the resulting metric turns out to be precisely the subcase 71 = of ([53l) below. 

If 7i 7^ 0, we can use Au/u as coordinate x, (as [5, A]z/ and [5, A]z/ show that 
{6 — (5)Az/ 7^ 0, z/, M and Au are certainly functionally independent). From (|46|1 we get 
X{u,x) = -l/(167iao), while ([43]) and ^ yield 

V 



U{u,iy,x) =- 

8ao7i 



( 2x^ - x%i + 2566a^7n + (52) 

Using a coordinate transformation x — > 16x71 (M)ao, and rescaling 71 and z/, we obtain 

1 



ds^ 
+2 



2n-kr^ -u + -il ( K^{h + x^fv"^ - 2A(6 - x^) + , ^ 



du 



dr + l^dz/ - 271 {K\h + x^y + 1) dx 



du + Az/W + ^. (53) 
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3.1.2. The class m{u) 
If m{u) ^ 0, we get, from ([45]) and ((471), that 

1 

dx 



sm?X 
m) and ([43]) then yield 

167iaz/ 



d'li dm 

Aau— h smuX — — 

du 



U 



X 



sm^ 



m[u) 

and ( 1461) can be rewritten as 

32 



ox^ ox ax 



Applying the coordinate transformation x 
transforms to 

dx J \ dy^ 

which gives 



y[x) 



(54) 

(55) 
(56) 

(57) 

with ^ = Xi, the latter equation 



+ s47] 



0, 



(58) 



(59) 



with S = sin or sinh according to the sign of s. We can now use the freedom in the 
choice of r to put Pi{u) equal to 1. 

After rescaling 71 and applying a transformation u mz//v^, we get the following 
expression for the line-element : 



ds^ 



2n - Ar^ - 4(s - Ai^'hf + 



+ 



2dr — 4z/vs — Az^7id?/ 



2m\/s — Ai/^ 
A 

2 971 



^7i 
dy 



du' 



Vs - Az/2 dy 



du 



du 



s- Az/2 



:du^ 



(60) 



3.2. The family n^O 
In this case ( l30l ) yields 



z/(^2 + A/3) 



22 



TT 



The Bianchi equations can be written as follows : 

(5v[/2 = -37rvlf2, 
= -37r\I^2, 
A^2 = 0, 
D^2 = 0, 



Sn 



7r((5z^ — utt) 

V 

-nibv — z/tt) 



Dtt 



V 
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while the NP Ricci equations yield 

3^2 + A - 67r2 



and 



'2 

a = , 7 = 7 

127r ' ' ' 



<5i. = _.''(2^2 + A) A7r = 0, 



;/(3^2 + A) 



Di^ = 0, D7 = -civj/f 



57 = 0, 
^7 = 0, 



We thus have the total derivatives : 

d7 = A-fiv^ - Cl^2^^a;^ 

z/(3^2 + A), 1 9x . ^ 

d\E'2 = -37r\l/2('*^^ + '^^)- 

Remember \1'2 = ^2(71") (l32ll . We can use vr (or \I'2) as a coordinate, but we prefer to 
write n = vr(x), \E'2 = ^2(2^) etc. 

The Cartan equations and the directional derivatives of tt show that we can write the 
basis one-forms as : 

u:^ = Sdx + i{Udu + Xdx + fdy), (61) 
LJ^ = du, (62) 
LJ^ = dr + Hdu + Fdx + Gdy, (63) 

where /, F, G, H, U and X are real functions of (m, r, x, y), and is a real function of x, 
which has to satisfy 
d^o 

= -657r^2, (64) 

or, equivalently, 

^ = -l5(67r2 + 3^2 + A). (65) 
dx 3 



The corresponding directional derivatives are given by 



^ ^iXG~fF + iSG)l: + f4--iX + iS)-^\, (66) 



25*/ \ dr dx dy 

D = ^. (68) 



Applying [A, D] and [5, D] to y shows /, U and X are independent of r, so a 
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transformation of y exists making [/ = 0. From [(5, A]?/ it then follows that / and 
X are independent of m so a further transformation of y exists, making X = 0. Next 
([32]) and ((61), show that 

/ = (69) 

a factor fof'{y) in / being absorbed in dy. From the directional derivatives of 7, 
we get 7 = — ci\l/2^^r + 'y'{u,x,y), F{u,r,x,y) = —4:7r{x)rS{x) + Fi{u,x,y) and 
G{u,r,x,y) = Gi{u,x,y), where 7' has to satisfy 

^ = -c^^rF^, (70) 

^ = -civl/f Gi. (71) 
dy 

Notice also that, in case Ci = 0, 7' is a function of u only, such that we can use the 
remaining boost freedom to put 7' equal to zero. This allows us to write 7' = cfb. 
The directional derivatives of z/, together with (l64l) . lead to an expression for u : 

n(u) , , 

nix)^l/\x) 



and from [A, D]r, [6, A]r and [6, 6]r we find an expression for H 
1 



H 



,(73) 



where m{u), n{u) and b{u,x,y) are real functions. 

If we use as coordinate x, and perform a translation of r — r + Cib(u, x,y) /x"^^^ , 

the line-element for this family is given by 

2 4r 
ds^ = — (cix'^/^r^ + nx^^^ — m) du^ + 2drdu H da;dM 

Notice that, in the case where Ci 7^ 0, we can fix the boost by requiring n = 1. 



4. Conclusions 



We can state that all Petrov type D pure radiation metrics, with or without cosmological 
constant, which admit a non-diverging null congruence, can be written in one of the three 
forms (l53l) . (I6OII or (j74l) . Metric (Tfil) reduces to the metric of Wils and Van den Bergh[3] 
for n{u) = 0. One can readily show that the subcase 71 = of (l53l) is a singular limit of 
(iGO]) . But whether or not (l53l) with 71 7^ is a singular limit of (l60|l is not yet clear. It is 
clear that in general ([741) admits at least one Killing vector field. More isometrics may 
arise for particular choices of the parameters and free functions. A detailed discussion 
of the physical properties of these metrics will be given in a follow up paper. 
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